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ABSTRACT 

Using photometric data from the 2MASS and GLIMPSE catalogues, I investigate 
the incidence of mid-infrared excesses {r^ 10 |J.m) of G and K stars of luminosity 
class III. In order to obtain a large sample size, stars are selected using a near-IR 
colour-magnitude diagram. Sources which are candidates for showing mid-IR excess 
are carefully examined and modelled to determined whether they are likely to be G/K 
giants. It is found that mid-IR excesses are present at a level of (1.8 ± 0.4) x 10^'^. 
While the origin of these excesses remains uncertain, it is plausible that they arise from 
debris discs around these stars. I note that the measured incidence is consistent with a 
scenario in which dust lifetimes in debris discs are determined by Poynting-Robertson 
drag rather than by collisions. 

Key words: infrared: stars - circumstellar matter. 



1 INTRODUCTION 

It has long been recognised (e.g. Jura 1999) that the presence 
of infrared excesses in giant stars which have yet to evolve 
, up the asymptotic giant branch may be due to the presence 
of long-lived debris discs around these stars. Recent work 
on debris discs has concentrated on their properties during 
I the main-sequence phase of stellar evolution but is becoming 
' clear that some debris discs survive into the post-main se- 
, quence lives of st ars, as exemplified by the whi t e dwarfs G29- 
■ 38 an d GD 362 (|Zuckerman fc Becklinlll987l : iBecklin et al.l 
' I2OO5I ) , a nd possibly ar ound the central star of the Helix 
Nebula (|Su et akl 120071 ). The presence of detectable discs 
around evolved stars may provide useful information about 
processes in debris discs since they have radiation environ- 
ments which are distinct from systems in which the host star 
is on the main-sequence. The study of debris discs might be 
especially fruitful for giant stars which have yet to ascend 
the asymptotic giant branch, since infrared emission is un- 
likely to be contaminated by emission from dust which has 
been lost from the star itself. 

To date, several studies have been made into the in- 
cidence of infrared excesses in giants. In the main, these 
studies have concentrated on the far-IR emission from gian t 
stars. An important result from iPlets fc Vvnckieij (|l999l ). 
which is based on det ailed analysis of da ta from the IRAS 
Faint Source Catalog (jMoshir et al.lll992l ). is that the inci- 
dence of 60 (i.m excesses in giants of spectral type G and K 
is 14 ± 5 per cent. 

In the mid-IR (i.e. wavelengths of order 10 H-m), the 
incidence of excesses in G/K giants is not well determined. 
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There have been indications that there are some such stars 
which do exhibit mid-IR excess e s. For example, a study by 
IClarke. Qudmaiier fc LumsdenI (|2005J), identifies four stars 
with spectral type G or K and with luminosity class III 
which show an excess at wavelengths of 8 (xm. Another 
example of a G/K giant which shows mid-IR excesses is 
HD 233517, which was subje ct to mid-IR spect roscopic ob- 
servations using Spitzer by Ijura et al.l l|2006l ). who con- 
cluded that this is likely to be a binary system with a pecu- 
har history. 

An indication that some G/K giants may have mid-IR 
excesses is also evident fr om a study of mid-IR extinction by 
llndebetouw et al.l (|2005l ). These authors noted, in passing, 
that in a selected sample of red-clump stars (which corre- 
spond to G/K giants) there appeared to be a small propor- 
tion of stars which exhibit excesses at 8 (xm in comparison 
to measured Ks magnitudes (see their figure 4). Indebetouw 
et al. attributed these mid-IR excesses to circumstellar dust 
but did not examine this population in detail. In particular, 
little attention was paid by these authors as to whether the 
stars with mid-IR excess are actually G/K giants. In this 
work I investigate a similarly selected population of stars, 
but take considerable care in the selection of such a pop- 
ulation and to pay particular attention to determining the 
nature of stars which show mid-IR excesses. 

The aim of this study is to determine the incidence of 
excesses at 8 \im in G and K giants, and to use available 
near- and mid-infrared photometry to investigate the nature 
of the sources showing such excesses. The most straightfor- 
ward approach to such a study would be to cross-correlate 
a large spectr al class catalogue, such as the Tycho-2 Spec- 
tral Catalog (|Wright et al.ll2003f) with a mid-IR cat alogue 
such as the GLIMPSE Catalog (JBeniamin et allbOOSh . This 
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approach was investigated but was found to be limited by 
the fact that the number of spectroscopically identified G/K 
giants returned by such a cross-matching is only of order of 
a few hundred, and of these, many are relatively bright and 
hence are saturated in one or more photometric bands of the 
Spitzer IRAC instrument. So, while this approach is satis- 
factory for investigating mid-IR exces ses in main-sequen ce 
stars (and is the approach adopted bv lUzpen et al.l l2007) it 
returns too few matches to provide any useful information 
on the incidence of such excesses in G/K giants. 

The approach adopted here is to use a near-IR colour- 
magnitude diagram to select a large sample (about 10*) of 
red-clump stars, and then to investigate in detail the stars 
from this sample which exhibit mid-IR excesses. In partic- 
ular, I attempt to ascertain whether their spectral energy 
distributions are consistent with being G/K giants with ex- 
cesses due to circumstellar dust and to rule out other plau- 
sible types of source. 

This paper is structured as follows; Section 2 describes 
the data sets and the survey region adopted for this study 
and how a sample of red-clump giants was selected. Section 
3 reports on the subset of these sources which are identified 
as having a mid-IR excess. The implications of the observed 
incidence of genuine mid-IR excesses in G/K giants is dis- 
cussed in Section 4. 



2 DATA SETS AND SURVEY REGION 

The basis of this study is to use a near-infrared colour- 
magnitude (CM) diagram to select a population of objects 
which have a high probability of being G or K giants. A dis- 
tinctive feature of any deep CM diagram of a sample which 
contains evolved stars of near-solar metallicity is the pres- 
ence of a feature arising from the so-called red-clump, which 
corresponds to the core-helium burning phase of intermedi- 
ate mass st ars. Studies of the propert ies of red-clump giants, 
such as by IZhao. Qiu fc Mad (|200l|) identify such stars as 
being of spectral type G or K and of luminosity class III (I 
shall use the terms red-clump stars and G/K giants inter- 
changeably here). Stars which can be identified as belonging 
to this population have been used in a variety of contexts as 
standard candles; for instance , to study Galactic structure 
l|L6pez-Corredoira et al.ll2002l ). In this study, the identifica- 
tion of red-clump stars is used as a basis for studying the 
properties of the stars themselves. In particular, the near-IR 
CM diagram is used to select a population of stars for which 
near and mid-IR photometry is available. 

A common choice of photometric bands which can be 
used to construct the required CM diagram is the near-IR J 
and K bands. In this study I use data f rom the 2MASS Point 
Source Catalog fSkrut skie et aLlbOOm . which comprises J, 
H and Ks photometric bands {Ks overlaps with the stan- 
dard K band, but is slightly narrower and has a centre-of- 
band wavelength of 2.15 nm). Figure 1 shows the (J — Ks) 
versus Ks colour magnitude diagram for the chosen survey 
region (a description of which is given below). In this dia- 
gram, the red-clump stars form the band of sources stretch- 
ing from (J — Ks) ~ 1.0, Ks ~ 11, down and to the right 
towards (J - Ks) « 1.8, Ks « 14. 

The identification of this feature with G/K giants is sup- 
ported by its proximity to the expected locus of red-clump 



stars given a reasonable model for extinction along lines of 
sight in the survey region. If it is assumed that the extinction 
per unit distance in the Ks band, ck, is uniform, then for 
a given value of ck , the tracks followed by different spectral 
cla sses of stars can be determined. Adopting the values given 
by iLopez-Corredoira et al.l l|2002r ) for canonical red-clump 
giants, that the absolute Ks magnitude Mk = 1-65 and the 
un-reddened colour {J — Ks)o = 0.75, it is found that a good 
match between the red-clump feature and the expected lo- 
cation of such stars is obtained for ck ~ 0.09 mag kpc~^ 
(indicated by the track labelled 'RCG' in Figure [T|. This 
va lue of Ck is in good a greement with the extinction data 
of [Marshall et al.l 1 20061 ) for the chosen survey region when 
determined over scales of several kiloparsecs. 

To the left and slightly below this feature (which I refer 
to as the red-clump band), the sources which have {J — Ks) 
in the range from to about 1, ar e likely to be main se- 
quence stars (|Bessell fc BrettlllQSg) . with later type stars 
showing relatively less reddening due to their proximity to 
the Sun. To the right of the red-clump band the CM dia- 
gram is populated by giants ascending the giant branch and 
the asymptotic giant branch. The degree of discrimination 
between red-clump giants and later type giants is illustrated 
in Figure [TJ which also shows the expected loci of stars of 
type MOIII and M5I II (with (J — Ks) colours as given by 
iBessell fc Brett 1988, and absolute magnitudes derived from 
values given by Cox 200 0) under the assumption of the same 
extinction model as adopted for the displayed locus of RCG 
stars. It can be concluded that the red-clump is a clearly 
defined feature in the colour-magnitude diagram. 

The survey area for this study was selected such that 
mid-IR p hotometry could be o btained from the GLIMPSE 
Catalog l|Beniamin et al.l l2003l fl. GLIMPSE is a survey of 
the Galactic plane (10° < |Z| ^ 65° and \b\ sj 1.0°) made 
using the four wavelength bands (3.6, 4.5, 5.8 and 8.0 yim) 
of the IRAC instrument on Spitzer l|Fazio et al.1 120041 '). A 
preliminary study for this work was based on the Version 1 
release of the GLIMPSE catalogues (released April 2005), 
while the results presented here are based on the Version 2 
release (April 2007) unless otherwise indicated. 

Any direction in the GLIMPSE survey region has a 
long line-of-sight through the disc of the Galaxy. This is 
advantageous in terms of the numbers of red-clump stars 
that can be detected, but, of course, the survey regions are 
subject to relatively high levels of extinction. Such extinc- 
tion, corresponding to column densities of order a few x 
lO'^^ H atoms cm~^, is not an issue as such, but it is im- 
portant to appreciate that variation of extinction across a 
survey region will lead to a loss of definition (i.e. a broad- 
ening) of the red-clump band in the resulting CM dia- 
gram. In choosing any extended survey area, some varia- 
tion of extinction is inevitable, but this variation can be re- 
duced by excluding regions with a high content of molecular 
gas. In choosing the survey region, the aim was to max- 
imise the area covered subject to the constraint of avoid- 
ing regions of high column density in molecular gas. This 
was done by avoiding known star formation regions, and 
by choosing a region of relatively low column density in 



^ All GLIMPSE data products were obtained from 
|http://www. astro. wisc.edu/sirtf/gllmpsedata.html| 
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Figure 1. The Ks versus [J — Ks) colour-magnitude diagram of the survey area. Sources included in this diagram are those which have 
photometric data in 2MASS at J, H, and Kg and GLIMPSE at 3.6, 4.5 and 5.8 ^un). The solid curves indicate the loci expected for 
red-clump giants (RCG) and M giants under the assumption of uniform extinction per unit distance in the Ks band of 0.09 mag kpc~^, 
with distances of 4 and 8 kpc indicated by dashed lines. The parallelogram (heavy lines) indicates the criteria used to select red-clump 
sources from this sample (see text for details). 



molecul ar hydroRen from inspection o f the CO emission 
map of iDame. Hartmann fc ThaddeusI (|200ll ) . The region 
used for this study is defined by 45.0° ^ (1) ^ 50.25° and 
0.70° ^ (b) ^ 1.0°. For the area covered by the GLIMPSE 
survey, the selected region has relatively low column densi- 
ties of CO. 

The number of 2MASS sources with measured values 
of J and Ks in the selected region is 93888. The extent of 
the survey area is such that about 10* clump giants would 
be included in the final sample, which is sufficient to mea- 
sure the incidence of mid-IR excesses at a level of order 0.1 
per cent. 



2.1 Selection of red-clump sources 

As noted above, the selection of red-clump sources is based 
on the (J — Ks)-Ks CM diagram. Rather than using all 
sources with 2MASS J and Ks data, preliminary studies 
revealed that an improvement in the degree of contrast be- 
tween the red-clump and the main sequence features can be 
achieved by requiring that included sources are detected the 
all 2MASS bands and the 3.6, 4.5 and 5.8 \uxi bands of the 
GLIMPSE Catalog. The number of sources in the survey 
region which satisfy this criterion is 53520 and it is these 
sources which are shown in the CM diagram in Figure [T] 

Stars which have a high probability of being G/K gi- 
ants were selected by following the leftmost side of the the 



red-clump band which is evident in Figure [T] While such 
a selection is somewhat subjective, using the left-hand side 
of the red-clump distribution in this way should minimise 
contamination of the sample from stars ascending the gi- 
ant branch/asymptotic giant branch. The selection criterion 
used was that Ks ^ 11.0 and that {J — Ks) should lie within 
±0.2 magnitudes of the straight line defined by 



(J - Ks) = 1.15 + 0.204(^-3 - 11.0) 



(1) 



This selection region is overlaid on Figure [T] Note that the 
cutoff at Ks = 11.0 was determined by the fact that at lower 
magnitudes the locus of the red-clump band becomes diffi- 
cult to trace due to the small number of sources in this part 
of the CM diagram. It can be seen from Figure [T] that this 
selection region is only an approximation to the expected 
track followed by red-clump stars under a model of uniform 
extinction per unit distance, but is a reasonable method of 
selecting such stars in the interval 11.0 < Ks < 13.5 (in 
fact, the faint cut-off adopted in this work is brighter than 
Ks ~ 13.5, see below). 

The CM diagram and selection region shown in Fig- 
ure[T]allows an estimate to be made of the maximum amount 
of contamination of the red-clump region by main- sequence 
stars. The profile of the distribution at constant Ks was 
examined at Ks ~ 13. The profile reveals a bimodal distri- 
bution with peaks due to main-sequence and evolved stars. 
At the minimum between the peaks of this distribution, the 
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density of main-sequence sources is about 10 per cent of 
the source density in the red-clump band. Since this mini- 
mum lies somewhat to the left of the red-clump band, and 
assuming that the density of main-sequence stars decreases 
as (J — Ks) increases, the contamination of main-sequence 
stars in the selected region will be less than 10 per cent: a 
reasonable estimate on an upper limit to the contamination 
would be at the level of a few percent. 

In addition to adopting selection criteria which isolate 
red-clump stars, it is also important to choose a faint cut-off 
in the value of Ks- This serves to reduce the bias inherent 
to detecting excesses at 8 \i.ra when many stars are close to 
the detection threshold at this wavelength. Ideally, it would 
be best to select a cut-off in Ks which results in all sources 
being detected at 8 vim. However, such a constraint leads 
to the rejection of a large proportion of the data. I investi- 
gated (using Version 1 of GLIMPSE) how the detectability 
of red-clump sources at 8 ti.m varies with Ks and found that 
at Ks = 12.9, 85 per cent of sources were detected. This 
cut-off in Ks was adopted as being a reasonable compro- 
mise between expected detection and source numbers. With 
these criteria applied (11 ^ Ks ^ 12.9 and Equation[T]) the 
number of red-clump sources selected is 9865. This sample 
forms the basis of the study of mid-IR excesses. 



2.2 The selection of sources with 8 ^mi excess 

An indication of the incidence of sources with excesses at 
8 \i.ra in the sample is revealed by the {Ks — [8.0])-(J — Ks) 
colour-colour diagram of the selected sample of red-clump 
sources (Figure [2^) . It can be seen that there is a small 
number (about 20) of sources with excesses in {Ks — [8.0]) 
which set them clearly apart from the main distribution of 
sources. Sources with less pronounced, but still significant 
excesses, can only be identified from more detailed analysis 
as described below. 

It can also be seen from Figure [2^ that the effects of 
interstellar reddening are apparent in the distribution of 
sources, with a trend for {Ks — [8.0]) to increase with in- 
creasing (J — Ks). It is therefore useful to define a colour 
{Ks — [8.0])o which includes the effects of reddening. The 
results of Indebetouw et al. (2005), and in particular, their 
best-fitting straight line to ([8.0] — Ks)-{J — Ks) data (as 
shown in their figure 4a and reproduced here on Figure [5^), 
allows this to be calculated as 

{Ks - [8.0])o = {Ks - [8.0]) - 0.37 x ( J - Ks) + 0.17 (2) 

For the sample of red-clump sources described above, the 
mean value, ({Ks-[S.O])o) is 0.048±0.002. The colour excess 
of any source is given by a colour excess defined by -Ea'8 ~ 
{Ks - [8.0])o - {{Ks - [8.0l)o) 

Sources which are candidates for stars having excesses 
at 8 p.m were selected using the criterion 

SNReks = -^ > 3.7 (3) 

The significance level here is determined by the sample size: 
it is chosen such that given a Gaussian distribution in Eks, 
it would be expected that there would be only one source 
in a sample size of 1 x 10"* which would exceed this value. 
On applying the above criterion, 52 sources were identified 
as candidate mid-IR excess sources. 



3 ANALYSIS OF CANDIDATE 8 nm EXCESS 
SOURCES 

In order to confirm or reject the candidate excess sources, 
the following approach was taken. Firstly, where doubts exist 
over the reliability of a source, it was culled from the sample. 
The remaining sources were then analysed in two ways: on 
the basis of their near-infrared colours and by modelling of 
the photometric data. 



3.1 Removal of low reliability sources 

The first stage in checking source reliability was to inspect 
the GLIMPSE Catalog source quality flag (SQF) for each 
source. Of the 52 candidate sources, 14 had SQF flags set 
which indicated confusion in in-band or cross-band merging. 
Two further sources were flagged as being potentially af- 
fected by banding correction and a single source was flagged 
as being possibly contaminated by stray light. 

In addition to the SQF data provided in the GLIMPSE 
Catalog, the 8 n.m GLIMPSE Image Atlas (Version 2) was 
visually inspected to check for artefacts which might have 
given rise to anomalous 8 ti.m fluxes. The highest resolu- 
tion (0.6 arcsec pixel size) 8 vxm GLIMPSE atlas images 
were used to check for anomalies. Fields of extent 1.2 ar- 
cmin around source positions were inspected and six sources 
(which would have been accepted on the basis of their SQF 
value of zero) were rejected on the basis of their appear- 
ance. Four of these sources appear to be affected by stray 
light contamination at 8 nm, while two are very close to rel- 
atively bright multiple sources. All six of these sources are 
also deemed to be of questionable reliability and, together 
with the 17 sources identified by their SQF values, were re- 
jected from the candidate list. The 2MASS and GLIMPSE 
photometry for the remaining 29 candidate sources is given 
in Table [T] Note that hereafter, individual sources are re- 
ferred to by the identifier (ID.) given in this table. 



3.2 Near-IR colours 

Some evidence as to the nature of the stars exhibiting 8 ^.m 
excesses can be determined by inspection of the (J — H)- 
{H — Ks) colour-colour diagram (Figure [^b) • The locations 
of the candidate 8 (im excess sources are indicated by the nu- 
merical identifier given in Table [l] The vast majority of the 
sources in this diagram are distributed in a band which cor- 
responds well to the expected location of G/K giants under 
different amounts of extinction (for reference, a reddening 
vector is shown which indicates the expected direction of 
the distribution of stars of type KOIII). While the majority 
of candidate 8 ti.m excess sources lie within the region popu- 
lated by reddened G/K giants, it is notable that seven of the 
candidate excess sources have (J — H) values which are sub- 
stantially lower than would be expected from their {H — Ks) 
colour (these are sources with identifiers 10, 13, 17, 18, 23, 
26 and 29). If these sources are G/K giants, they would ap- 
pear to have anomalous near-IR colours. The remaining 22 
candidate sources appear to be consistent with being G/K 
giants. Further investigation (below) uses the available 7- 
band photometry available to model these sources. 
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Figure 2. Colour-colour diagrams(a) {Kg — [8.0])-(J — ii's), and (b) (J — H)-(H — Kg), for all of the selected red-clump stars (grey 
asterisks). Candidate sources with excesses at 8 fim are indicated by the numerical identifier as given in Table [T] The straight line in (a) 
shows the extinction relation as derived bv llndebetouw et al.l BOOST) using a similar sample of stars. The reddening vector shown in (b) 
starts at the expected position of an un-reddened star of type KOIII. 
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3.3 Modelling of 7-band photometry 

The 29 sources with excesses in their (Ks — [8.0])o colours 
were further analysed to attempt to determine if these 
sources were likely to be G/K giants with excesses at 8 \ira, 
or if they have some other origin. As noted above, the sample 
of red-clump stars is expected to be contaminated by other 
types of source at the level of maybe as much as a few per- 
cent. However, it could plausibly be the case that it is these 
contaminating stars which give rise to the sources detected 
with 8 \ira excesses. Hence it is important then to consider 
whether the detected excess sources could be other types 
of star. The approach taken was to use the data from the 
seven photometric bands (J, H, Ks, 3.6^m, 4.5 i^m, 5.8 \im, 
and 8.0 \im) to compare to the spectra from a selection of 
plausible models. 

Clearly, one model which needs to be tested is that of a 
red-clump giant with a mid-IR excess. As a reference point, 
it is useful to define a 'ca nonical' red - clump giant using 
the results of the study of IZhao et al.l l|2001f ). The subset 
of metal rich ([Fe/H] > —0.3) clump giants in their study 
have a mean effective temperature of Tgg- = 4860 K and 
surface gravity of logg = 2.81. The canonical red-clump gi- 
ant can be approximated by the closest Kurucz model (see 
below), which has T^g — 4750 K, and logg — 3.0. 

In model fitting however, I adopted a range of giant star 
models with a range of surface temperatures from 3500 K 
to 6500 K. This is a wider temperature range than would be 
expected for red-clump giants, but provides a useful check 
that individual stars are not better modelled as cooler or 
hotter stars. This, and all other photospheric emission (in 
other models des cribed below) was model led using a Kurucz 
spectral modeO (|Cci8telli fc Kurucal2004r ) with solar metal- 
licity (also, the particular set of models used have a tur- 
bulent velocity of 2 km s^^). The temperature resolution of 
the published grid of Kurucz models is 250 K. 

The mid-IR excess is modelled as a component which 
only contributes significantly to the 8 yun band. This was 
implemented as a Gaussian component at a wavelength of 
7.9 nm with a — 0.6 ^xm, with only the normalisation of 
this parameter being allowed to vary. While this component 
is a reasonable representatio n of a PAH featur e observed in 
the spectrum of HD 233517 (|jura et all 120061 ') it should be 
stressed that an excess in one photometric band could be 
modelled in a variety of ways, and that I do not claim any 
particular significance to modelling the excess in this way. 

A quite difi'erent set of stars which are known to show 
mid-IR excesses are hot stars with infrared excesses aris- 
ing from bre msstrahlung emission from pla sma surround- 
ing the star IjGehrz. Hackwell fc Joneall974). Such sources 
are detected in considerable numbers bv lClarke et al.l (|2005l ) 
and Uzpen et al. (2007) and are believed to occur for stars 
of spectral type B8 or earlier. Consequently I also examine 
a hot star model comprising a main-sequence B8 star and 
bremsstrahlung emission from a plasma at a temperature of 
1 X 10*K. Note that as far as fitting to the combined 2MASS 
and GLIMPSE photometry is concerned, the model is rather 
insensitive to the exact stellar type: the choice of B8V is a 
conservative one which represents the coolest and lowest lu- 



minosity star that might plausibly also exhibit an infrared 
excess due to bremsstrahlung emission. 

A component of all of the fitted models was the in- 
clusion of extinction due to interstellar dust. The extinc- 
tion values used were taken from the publicly available 
synthetic extinction curvejj which are based on the grain 
m odel of iLi fc Draind (12001 ) an d the grain size distribution 
of IWeingartner fc Draind 1 20011 ). The particular extinction 
curve used was the one for which Rv = 3.1 and which these 
authors consider to provide a good match to the extinction 
properties of diffuse HI clouds in the Milky Way. 

For all of the models, three parameters were fitted: the 
normalisation X of the stellar spectral component, the nor- 
malisation of either the bremsstrahlung component or of the 
Gaussian line at 7.9 ^m, and the line-of-sight column density 
Nh- In addition to these parameters, note that for the red- 
clump giant model, the fitting was conducted across a range 
temperatures (as described above) in order to also deter- 
mine the photospheric temperature. A spectrophotometric 
distance d to the star can be derived from the normalisa- 
tion parameter X hy d = rtX~^'^ where r, is the stellar 
radius. Stellar radii were estimated by assuming blackbody 
emission and relationsh ips between stellar luminos i ty an d 
temperature as given bv lde Jager fc NieuwenhuiizenI l|l987h . 

It is important to appreciate that the spectral resolu- 
tion provided by the seven photometric bands used here is, 
in general, insufficient to provide any clear discrimination 
in luminosity class: data which are well modelled as a G/K 
giant with a Gaussian excess at 8 ^.m, will also be accept- 
ably fitted by models in which the star is a G/K dwarf or 
supergiant. Models in which the star is modelled as a dwarf 
and as a supergiant were also fitted to obtain spectropho- 
tometric distances and values of column density for the ex- 
cess sources. In order to investigate whether these models 
could be plausible, a set of neighbouring (non-excess) red- 
clump giants was used to estimate the column density in 
the direction of each candidate excess source. The set of 
neighbouring sources for each candidate excess source was 
defined as those sources which were selected according to the 
red-clump criterion (Equation [ij and which lie within 5 ar- 
cmin of the position of the excess source and which have Ks 
within 0.2 mag of this source. In fitting to these neighbour- 
ing red-clump giants, the canonical red-clump stellar model 



was assumed (i.e. T, 



eff 



4750 K, and logg = 3.0). The 



fitting process returned a column density and flux normali- 
sation {X) for each source, allowing a mean column density 
and mean spectrophotometric distance for the neighbour- 
ing red-clump stars to be determined. These values for a 
so-called 'neighbouring field' are given in Table [5] and are 
discussed further in Section 3.4. 

For a particular model, the expected flux density in each 
photometric band was obtained by integrating uniformly 
over a passband between the wavelengths at which the fll- 
ter response is half its maximum value. Such an integration 
is an approximation to modelling the exact photometric re- 
sponses of the filters, but is expected to be adequate for this 
work. 

The fitting process was a minimisation of chi-squared 
type approach, using the MPCURVEFIT routine in IDL. 



Obtained from|http://wwwuser. oat. ts. astro. it/castelli/grids.html| ^ http://www. astro. princeton.edu/~'draine/dust/dustmix. html 
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Best fit model 






Neighbouring field 


ID. 


Designation 


Type 


T/K 


xl 


Nh 


d/kpc 


Es 


SNRg 


Nh 


d/kpc 


1 


G045. 1000+00.8808 


RCG 


5750 


0.69 


12.7 + 0.5 


2.1 + 0.4 


0.13 


2.3 


8.3 + 1.0 


4.4 


2 


G045. 1844+00.7612 


RCG 


5000 


0.33 


9.8 + 0.9 


4.4 + 1.0 


0.89 


9.2 


9.1 + 1.1 


5.4 


3 


G045. 2769+00.7123 


RCG 


4500 


0.40 


7.9 + 0.9 


6.7 + 1.7 


0.52 


3.9 


9.1 + 1.1 


5.4 


4 


G045. 2852+00.9988 


RCG 


4500 


0.41 


7.7 + 1.0 


7.0 + 2.1 


0.52 


5.0 


9.2 + 0.8 


5.5 


5 


G045. 4108+00.8768 


RCG 


4500 


0.38 


5.8 + 0.9 


5.9 + 1.5 


0.42 


5.7 


7.4 + 2.4 


4.6 


7 


G046. 2397+00.9526 


RCG 


5000 


0.13 


6.6 + 0.6 


2.7 + 0.4 


0.28 


5.9 


6.8 + 1.3 


3.3 


8 


G046. 3252+00.8710 


RCG 


4750 


0.58 


9.6 + 0.7 


5.4+1.1 


0.33 


4.2 


9.2 + 1.1 


5.5 


9 


G046. 6118+00.9688 


RCG 


4500 


0.81 


8.0 + 0.9 


5.2 + 1.3 


0.27 


4.0 


8.3 + 0.9 


4.1 


10 


G046. 6781+00.7324 


HSB 


11500 


0.66 


15.2 


> 1.8 


2.45 


82 






11 


G046. 6833+00.8797 


RCG 


4500 


0.90 


9.3 + 0.9 


6.1 + 1.5 


0.35 


4.6 


9.6 + 0.9 


4.9 


12 


G046. 7233+00.7512 


RCG 


4250 


1.05 


7.9 + 0.8 


8.6 + 2.2 


0.54 


4.5 


9.8 + 0.7 


5.2 


13 


G046. 9192+00.7417 


HSB 


11500 


0.65 


13.1 


> 1.1 


0.65 


19 






14 


G047.5359+00.8104 


RCG 


5750 


1.96 


8.3 + 0.6 


1.8 + 0.3 


0.12 


2.7 


7.6 + 1.0 


3.3 


15 


G047.8477+00.9287 


RCG 


4000 


1.36 


4.9 + 0.6 


8.6 + 3.2 


0.64 


5.4 


7.6 + 1.0 


4.0 


16 


G048. 1801+00.9573 


RCG 


4250 


0.26 


6.9 + 1.1 


8.2 + 2.2 


0.43 


4.3 


8.4 + 0.9 


5.1 


17 


G048. 3706+00.7076 


HSB 


11500 


1.24 


15.3 


> 1.7 


0.75 


11 






19 


G048. 3924+00.8103 


RCG 


5000 


0.78 


12.1 + 1.1 


4.4+1.5 


1.04 


12 


9.5 + 1.3 


5.5 


20 


G048. 4911+00.8260 


RCG 


4500 


0.96 


9.2 + 1.1 


6.2 + 1.9 


0.57 


4.0 


9.0+1.1 


5.1 


21 


G048. 5964+00.7855 


RCG 


4500 


0.15 


9.2 + 0.8 


6.5 + 1.6 


0.40 


3.9 


9.7+1.1 


5.2 


22 


G048. 5970+00.9849 


RCG 


4750 


0.74 


7.2 + 0.8 


5.5 + 1.2 


0.40 


4.0 


8.4 + 1.1 


5.2 


23 


G048. 6036+00.9284 


HSB 


11500 


0.51 


17.7 


> 1.6 


0.75 


12 






24 


G049. 0793+00.7610 


RCG 


4500 


0.57 


5.0 + 1.0 


5.6 + 1.6 


0.33 


4.6 


7.4+ 1.2 


4.3 


25 


G049. 1099+00.8232 


RCG 


4500 


0.59 


7.0 + 0.9 


6.5 + 1.6 


0.48 


5.2 


8.7+1.1 


5.1 


26 


G049. 4210+00.8476 


HSB 


11500 


0.60 


13.5 


> 1.6 


0.76 


11 






27 


G049. 5812+00.9666 


RCG 


4250 


0.30 


6.8 + 1.0 


8.6 + 2.9 


1.23 


7.2 


8.0 + 0.9 


5.4 


28 


G049. 9290+00.8456 


RCG 


6250 


0.32 


14.0 + 0.6 


1.8 + 0.3 


0.23 


3.0 


9.2 + 1.0 


5.6 


29 


G050. 1256+00.9808 


HSB 


11500 


0.57 


16.3 


> 1.6 


1.15 


26 







Table 2. The results of model fitting. Column 1 shows the source identifier (ID.) as used in Table[T] Column densities (Nh) are expressed 
in units of (lO^^H atoms cm~^). The quantities Eg and SNRg are defined in the text. 



Given that fitting was carried out at fixed photospheric tem- 
peratures, the uncertainties in the fitted parameters which 
were returned by this process were unrealisticalfy small. 
Consequently, uncertainties in fitted parameters (and quan- 
tities derived from these parameters) were estimated by in- 
spection of the parameters returned from a range of temper- 
atures corresponding to variation in x^ oi about 1. 



3.4 Results of modelling 

The modelling process has two distinct aims. The first is to 
determine whether a given excess source can be acceptably 
fitted by a model which comprises a red-clump giant and 
Gaussian excess at 8 (xm (a model referred to here as RCG) 
or by a model comprising a hot star and bremsstrahlung 
emission (which is referred to as HSB). The second aim is to 
determine which of the those stars which are well modelled 
by the RCG model do actually exhibit a significant excess 
at 8 [i-m. 

Of the 29 candidate sources, 27 could be fitted by either 
the RCG or HSB models within the 95 per cent confidence 
interval and are taken to be acceptably modelled. The two 
sources (6 and 18) which were not acceptably fitted had min- 
imum values of x2 > 5 and hence the models can be rejected 
at the 99.9 percent confidence level. The spectral energy dis- 
tribution of source 6 (shown in Figure [S] left-hand panel) 
suggests that it has a continuum level which differs between 
the 2MASS and GLIMPSE observations. Since 2MASS and 
GLIMPSE data were collected at different times, a likely ex- 
planation is that this is a variable star and it is discounted 



from any further analysis. The nature of source 18 is less 
clear: this could be an RCG source with infrared excesses at 
as short a wavelength as 4.5 (xm, although its spectral en- 
ergy distribution (Figure [S] right-hand panel) also appears 
to be reasonably well-described by the HSB model with the 
exception of the flux measurement at 3.6 ^im. While it is 
possible that variants to the RCG or HSB model could de- 
scribe the data for this source, such modelling is not justified 
using these data. It is also noted that this source is one of the 
seven sources in Figure [2]d which lie away from the expected 
distribution of G/K giants (see Section 3.2). The nature of 
source 18 is acknowledged to be uncertain and this source is 
excluded from further analysis. 

For all of the sources which can be acceptably fitted 
by either the RCG or the HSB models, an excess at 8 (xm 
(called Es) was determined by calculating the difference be- 
tween the observed value of [8.0] and the model photospheric 
value ([8.0]phot)- The signal-to-noise ratio of the excess was 
calculated as 



SNRg = 



[8.0] - [8.0]phot 



'■[8.0] 



(4) 



This quantity is tabulated, along with the results of the 
model fitting in Table [2] 

Of the 27 sources which are well modelled by either 
the HSB or RCG models, it should be noted that the dis- 
tinction between these two classes of models is very strong: 
no source can be acceptably modelled by both types of 
model. Six of the 27 sources are better fitted by the HSB 
model than by the RCG model, and all of these sources 
were identified from Figure [2] as being unlikely to be G/K 
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SSTGLMC_045. 6806 + 00. 7488 



SSTGLMC_048. 3865 + 00. 7629 




Wavelength / /it 




Wavelength / /it 



Figure 3. The spectral energy distributions of the two candidate excess sources (6 and 18, shown on left and right respectively) which 
are not well fitted by either the RCG or HSB models (see text). Photometric data points are indicated by 'x' with error bars. For 
comparison the spectrum of a canonical red-clump giant is also shown (continuous line) together with the model photometric points 
(grey bars - the width of which represent the FWHM of the photometric bands). 



giants. The sources are characterised by high column den- 
sities (> 10^'^H atoms cm~^) and very strong mid-IR ex- 
cesses (SNRs > 10). Note that because the model assumes 
the star is the lowest luminosity star which is likely to ex- 
hibit bremsstrahlung emission (i.e. type B8V) the distances 
quoted in Table [2] are lower limits. Two examples of the 
spectral energy distribution of sources that are best fitted 
by the HSB model are shown in Figure|4] It is noticeable that 
this model can reproduce a wide range of spectral slopes at 
A > 2 fi.m. In particular, the source 10 seems to have spec- 
tral energy distribution which rises between the 2MASS and 
the GLIMPSE bands, and it is clear that this behaviour can 
be successfully described by the HSB model. 

The remaining 21 sources are well-fitted by the RCG 
model, and it is notable that in Figure O all of these sources 
are at locations consistent with them being G/K giants. The 
quantity SNRg is used to determine whether the excess at 
8 (i.m is confirmed. As in the case for selecting candidate 
excess sources using Eks,, the cut-off is at a level where only 
one excess source is expected to occur by chance from the 
total sample of red-clump giants, i.e. SNRg > 3.7. Of the 
21 candidate excess sources which can be fitted by the RCG 
model, 18 show excesses which satisfy this criterion, and of 
these, seven sources have SNRg ^ 5. 

Three sources (1, 14, and 28) do not have confirmed 
excesses. It is noticeable that these are all of high photo- 
spheric temperature (Tgfj- = 5750 K, 6250 K). Such temper- 
atures would be unusually high for a red-clump star, and it 
is possible that these stars are low metallicity stars on the 
horizontal branch. However, given that their mid-IR excesses 
are unconfirmed, they are rejected from the sample. 

The remaining 18 sources are strong candidates for be- 
ing G/K giants with mid-IR excess. Examples of the spec- 
tral energy distributions of two of these sources are shown 
in Figure [5] As might be expected from the goodness-of- 
fit for these sources, the RCG model provides an excellent 
match to these data, with no indication of systematic devi- 
ations which might cause us to question the validity of the 
model. However, one possibility that should be addressed is 



that these stars might be of a different luminosity class to 
that assumed in the model, i.e. that the stars may be on the 
main-sequence or supergiants. 

This was investigated by fitting using Kurucz models 
corresponding to G/K dwarfs and supergiants (using Ver- 
sion 1 of the GLIMPSE data) . In the case of models which 
are well-fitted as red-clump giants, the best-fitting dwarf 
or supergiant model typically has a very similar tempera- 
ture (usually differing by no more than 250 K) and similar 
column density to the best-fitting giant star model. The dis- 
tance estimates, as would be expected, depend very strongly 
on the assumed luminosity class. As can be seen from Table[2] 
the distances implied by giant star models are broadly con- 
sistent with those obtained from the neighbouring red-clump 
stars (labelled as the 'neighbouring field' in Table [5)|. Dwarf 
star models (luminosity class V) imply distances which are 
typically an order of magnitude lower than the mean dis- 
tance to the neighbouring red-clump giants, but at a similar 
value of total column density. It is concluded that it is highly 
unlikely that any of the candidate excess stars are on the 
main-sequence. The distance estimates which are arrived at 
assuming that the candidate excess star is a supergiant of 
luminosity class I are typically of order 100 kpc. The as- 
sumption that these stars are supergiants can be rejected 
since such distances would place these stars well outside the 
stellar disc of the Milky Way. It seems reasonable to conclude 
that the candidate excess stars are all giants (i.e. luminosity 
class HI). 



4 DISCUSSION 

The analysis presented here confirms that there is a pop- 
ulation of red-clump stars which appear to be G/K giants 
with significant mid-IR excesses. Out of an initial sample 
of 9865 red-clump stars, the analysis presented here identi- 
fies 18 sources which fall into this category. The incidence 
of mid-IR excesses in G/K giants is thus calculated to be 
(1.8 ± 0.4) X 10~^. This is substantially lower than the in- 
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Figure 4. The spectral energy distributions of two excess sources (10 and 13, sliown on left and right respectively) which are well fitted 
by HSB models (see text). These two examples illustrate the wide range of continuum slope which can be generated by this model. The 
spectral energy distributions are displayed in the same manner as in Figure [3] 
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Figure 5. The spectral energy distributions of two excess sources (2 and 27, shown on left and right respectively) which are well fitted 
by RCG models (see text). These are the two sources of this type with the strongest 8 \xm excesses. The spectral energy distributions 
are displayed in the same manner as in Figure |3] 



cidence of f ar-IR excesses in G/K g iants of 14 ± 5 percent 
reported bv lPlets fc Vvnckierl (|l999r ). This low-incidence of 
mid-IR excesses as compared to far-IR excesses mirrors, in 
a qualitative sense at least, the behaviour seen for main- 
sequence stars, where mid-IR excesses are far less com- 
mon than those detected at far-IR wavelengths, which is 
interpreted as reflecting the incidence of systems with hot 

i '^ few X 100 K) and cold (< 100 K) dust respectively 
Uzpen et al.ll2007l ). A discussion of implications of our ob- 
served incidence of mid-IR excesses in G/K giants is given 
below. 

In investigating candidate excess sources, this work has 
also identified a substantial number of sources whose SEDs 
are well-described as highly absorbed early-type stars in 
which the excess emission arises from bremsstrahlung emis- 
sion (the HSB model). With the admittedly limited infor- 
mation provided by the 7-band photometry used here, some 
caution should be expressed about the physical interpreta- 



tion of this model: it may plausibly be the case that this 
type of SED is produced by a G/K giant (or other luminous 
star) which has infrared excesses across most of the photo- 
metric bands used here. An argument against this view is 
however provided by the (J — H)-{H — Ks) diagram (Fig- 
ure (2]). When the suspected variable star (source 6) and the 
high-temperature RCG sources (1, 28 and 14) are removed 
from the diagram, the distinction between RCG and HSB 
models becomes quite clear, with the HSB stars seeming 
to follow an extinction vector which is distinct from that 
followed by G/K giants. If the stars which are here classi- 
fied as HSB were simply RCG stars with different types of 
infrared excess, then such a distinct separation would not 
be expected. So it seems likely that the stars modelled as 
HSB type are distinct from G/K giants, although without 
spectroscopic identification of the stellar types, it cannot be 
confirmed that the hot-star plus bremsstrahlung model pro- 
vides the correct physical interpretation for these sources. 
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It is also interest to examine whether the number of 
HSB sources detected here is consistent with the known 
frequency of early-type stars with bremsstrahlung emis- 
sion. This work reveals six such sources in a survey area 
of 1.575 deg^, i.e., at a surface density of about 4 deg~^. 
Source s of this typ e are ubiquitous in studies of mid-IR 
excess: lUzpen et al.l (120071) report 126 such sources, which 
corresponds to a surface density (in the GLIMPSE survey 
region) of about 0.6 deg~^. In that work however, the typ- 
ical Ks magnitude of the early-type star involved is about 
8. The work presented here is sensitive to Ks ^ 12.9. A 
crude extrapolation (based on a uniform spatial distribu- 
tion of sour ces) of the su r face d ensity of these sources from 
the work of lUzpen et al.l (|2007l l to the flux limits adopted 
here yields a value of order 100 sources deg~ . It does not 
seem unreasonable then that a small fraction of such sources 
fall within the region of the Ks/{J — Ks) colour magnitude 
diagram used to isolate RCG stars. 

Turning now to the sources which are identified as be- 
ing G/K giants with excesses at 8 ^im, a key question relates 
to the nature of the observed mid-IR excess. If the emission 
arises from a debris disc, it might be expected that the spec- 
trum of the emission would be reasonably well approximated 
by a blackbody spectrum. Given that the photometric data 
presented here do not extend to wavelengths beyond 8 ^i.m, 
it is not possible to determine a blackbody temperature of 
the component giving rise to the mid-IR excess. However, 
an estimate can be made on the upper limit to the black 
body temperature, by requiring that the excess component 
does not exceed measured values at 5.8 ixm. Inspection of 
blackbody models in comparison to the SEDs of two strong 
excess sources (2 and 27, whose SEDs are shown in Fig- 
ure [5]) reveals that models in which mid-IR emission arises 
from a blackbody component are a good representation of 
the data when the blackbody temperatures are lower than 
about 450 K, and temperatures in excess of 600 K can clearly 
be ruled out. 

Another possible source of mid-IR emission is from 
PAHs. This is particularly pertinent to this study as one of 
the few G/K giants w hich has been sub ject to mid-IR spec- 
troscopy, HD 233517 (|jura et al.ll2006h . shows strong PAH 
emission. In the 8 ji-m band, the infrared excess of this K2III 
star is dominated by PAH emission which has an ampli- 
tude which is approximately 70 percent of the photospheric 
emission in this band (estimated from figure 1 of Jura et 
al.). The model used here to describe red-clump giants with 
mid-IR excesses included a Gaussian component which was 
chosen to provide a reasonably good representation of the 
7.9 p.m PAH feature in the spectrum of HD 233517. I note 
that several sources have 8 ti.m excesses of about 0.5 magni- 
tudes, which corresponds to the relative amplitude of Gaus- 
sian component with respect to the photospheric emission 
being similar to that observed in HD 233517. On the basis 
of these data I would not claim that PAH emission is the 
dominant cause of mid-IR excesses all G/K giants, but note 
that it is a possibility which deserves further investigation. 

It is also of interest to examine the consequences of the 
measured incidence of mid-IR excesses in G/K giants in re- 
lation to the corresponding incidence in their main-sequence 
precursors. It is known that main-sequence stars give rise to 
mid-IR excesses which most likely arise from hot dust in de- 
bris discs. Dust is produced by collisions within debris discs 



but its emission depends on the radiation field of the host 
star. As a star with a disc evolves beyond the main-sequence, 
it is likely that the collisional processes are relatively unaf- 
fected, whereas the radiation environment changes as the lu- 
minosity (and temperature) of the star change. An increase 
in luminosity of the host star between the main-sequence 
and the red-clump has two distinct effects. The first is that 
the blackbody temperature of dust (at a given radius from 
the star) will increase by a factor (I/rc/-I'Ms)^ , where Lms 
and Lrc are the stellar luminosities on the main-sequence 
and in the red-clump respectively. Secondly, if the lifetime 
(tpr) of dust (i.e. from its generation to its accretion onto, or 
sublimation close to, the star) is determined by Poynting- 
Robertson (P-R) drag, then this lifetime will vary according 
to 



ipr.RC _ Lms 

ipr,MS -^RC 



(5) 



(where the subscripts MS and RC refer to main-sequence 
and red-clump respectively). While the assumption of the 
importance of P-R drag would be contrary to the theoretical 
expectation that dust lifetimes are deterrnined b y collisional 
timescales in massive debris discs l|Wvattll2005l ). it is useful 
to explore such a scenario in the light of the results presented 
in this paper. 

In order to analyse how these changes might have an 
impact on the observations presented here, it is necessary to 
identify the main-sequence precursors to red-clump stars. 
While recogiusing that red-clump stars, can in principle, 
have a relati vely wide r ange of masses (O.5M0 to 2.5M0, 
see e.g. IZhao et al.ll200H ) I choose a value of 2Mq as being 
representative of the mass of the main-sequence precursor. 
This would correspond to spectral type A5V, with a lumi- 
nosity of about 121/0 . If such a star evolves to form a canon- 
ical red-clump star, its luminosity would increase to about 
46L0. 

The effect that this would have on dust temperature is 
that it would increase by about 40 per cent. The detection 
of a difference in mean temperatures between a population 
of main-sequence stars and red-clump stars would require 
a much more comprehensive survey of mid-IR spectra than 
are currently available and certainly would require longer 
wavelength data than used in this study. Even if such a 
survey could be conducted, the known variation in main- 
sequence disc temperatures, ranging from about 200 K to 
800 K (Uzpen et al., 2007) would make it very difficult to 
interpret such a study. 

On the other hand, the effect of variation of stellar lu- 
minosity on fpr may have more readily apparent observable 
consequences. Under the assumption that dust lifetimes are 
dominated by P-R drag, then provided that dust is gen- 
erated by episodi c events (such as in the solar system, e.g. 
iFarlev et al.ll2006r ). and that the mean interval between dust 
generating events is much greater than ipr, then the inci- 
dence of mid-IR excesses will be proportional to tpr. 

The observed incidence of mid -IR excesses i n A t ype 
main-sequence stars is reported by I Uzpen et al.l l|2007h as 
1.0 ± 0.5 per cent. Using the values given above for Lms 
and I/Rc, the ratio of the decay timescales tpr.Rc/ipr.MS = 
0.27, and within the scenario described above, this should 
correspond to the ratio of incidences of mid-IR excesses. 
Hence the expected incidence of mid-IR excesses in G/K 
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giants would be expected to be (2.7 ±1.4) x 10"^. This is 
in agreement with our observed value of (1.8 ± 0.4) x 10"'^. 
Hence, a model in which dust lifetimes in debris discs are 
dominated by P-R drag leads to a natural explanation of 
the observed incidence of mid-IR excesses in G/K giants. 



5 CONCLUSIONS 

The major conclusion of this study is that mid-IR excesses 
are exhibited by G/K giants albeit at a very low level of 
incidence. This could be consistent with the level of mid- 
IR excesses in their main-sequence precursors provided that 
dust generation is episodic and dust lifetimes are determined 
by P-R drag rather than collisional processes. In passing, I 
note that this work has obtained a value for the incidence of 
mid-IR excesses in G/K giants which is better determined 
than the corresponding incidence in main-sequence stars. 
While it is assumed here that the mid-IR emission arises 
from relatively large dust grains in a debris disc, it should 
be emphasised that the nature of the mid-IR emission in this 
population remains uncertain, and requires mid-IR spectro- 
scopic investigation. 

This study has illustrated that the technique of isolat- 
ing samples of red-clump stars from CM diagrams can be a 
useful technique to study this population of stars. I would 
expect that this technique could be put to good use by its 
application to a study of the incidence of excesses at longer 
infrared wavelengths. 
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Table 1. Photometric data for reliable sources selected as candidates for having 8 ^rn excesses. Column 1 (ID.) is the identifier used to refer to individual sources in this paper. For 
each photometric band a() represent the single standard deviation uncertainties. The quantities Eks and SNRbk8 ^I'e defined in the text. 



